Multimodal MRI study on the relation between WM integrity and connected GM atrophy and its effect on disability in early multiple sclerosis.

BACKGROUND: Multiple sclerosis (MS) is characterized by pathology in white matter (WM) and atrophy of grey matter (GM), but it remains unclear how these processes are related, or how they influence clinical progression. OBJECTIVE: To study the spatial and temporal relationship between GM atrophy and damage in connected WM in relapsing-remitting (RR) MS in relation to clinical progression. METHODS: Healthy control (HC) and early RRMS subjects visited our center twice with a 1-year interval for MRI and clinical examinations, including the Expanded Disability Status Scale (EDSS) and Multiple Sclerosis Functional Composite (MSFC) scores. RRMS subjects were categorized as MSFC decliners or non-decliners based on ΔMSFC over time. Ten deep (D)GM and 62 cortical (C) GM structures were segmented and probabilistic tractography was performed to identify the connected WM. WM integrity was determined per tract with, amongst others, fractional anisotropy (FA), mean diffusivity (MD), neurite density index (NDI), and myelin water fraction (MWF). Linear mixed models (LMMs) were used to investigate GM and WM differences between HC and RRMS, and between MSFC decliners and non-decliners. LMM was also used to test associations between baseline WM z-scores and changes in connected GM z-scores, and between baseline GM z-scores and changes in connected WM z-scores, in HC/RRMS subjects and in MSFC decliners/non-decliners. RESULTS: We included 13 HCs and 31 RRMS subjects with an average disease duration of 3.5 years and a median EDSS of 3.0. Fifteen RRMS subjects showed declining MSFC scores over time, and they showed higher atrophy rates and greater WM integrity loss compared to non-decliners. Lower baseline WM integrity was associated with increased CGM atrophy over time in RRMS, but not in HC subjects. This effect was only seen in MSFC decliners, especially when an extended WM z-score was used, which included FA, MD, NDI and MWF. Baseline GM measures were not significantly related to WM integrity changes over time in any of the groups. DISCUSSION: Lower baseline WM integrity was related to more cortical atrophy in RRMS subjects that showed clinical progression over a 1-year follow-up, while baseline GM did not affect WM integrity changes over time. WM damage, therefore, seems to drive atrophy more than conversely.

Effects of Acquisition Parameter Modifications and Field Strength on the Reproducibility of Brain Perfusion Measurements Using Arterial Spin-Labeling.

BACKGROUND AND PURPOSE: Although the added diagnostic value of arterial spin-labeling is shown in various cerebral pathologies, its use in clinical practice is limited. To encourage clinical adoption of ASL, we investigated the reproducibility of CBF measurements and the effects of variations in acquisition parameters compared to the recommended ASL implementation. MATERIALS AND METHODS: Thirty-four volunteers (mean age, 57.8 ± 17.0 years; range, 22-80 years) underwent two separate sessions (1.5T and 3T scanners from a single vendor) using a 15-channel head coil. Both sessions contained repeated 3D and 2D pseudocontinuous arterial spin-labeling scans using vendor-recommended acquisition parameters (recommendation paper-based), followed by three 3D pseudocontinuous arterial spin-labeling scans, two with postlabeling delays of 1600 and 2000 ms and one with increased spatial resolution. All scans were single postlabeling delay. Intrasession (identical acquisitions, scanned five minutes apart) and intersession (first 2D and 3D acquisitions of two sessions) reproducibility was examined as well as the effect of parameter variations on CBF. RESULTS: Intrasession CBF reproducibility was similar across image readouts and field strengths (within-subject coefficient of variation between 4.0% and 6.7%). Intersession within-subject coefficient of variation ranged from 6.6% to 14.8%. At 3T, the 3D acquisition with a higher spatial resolution resulted in less mixing of GM and WM signal, thus decreasing the bias in GM CBF between the 2D and 3D acquisitions (ΔCBF = 2.49 mL/100g/min [P < .001]). Postlabeling delay variations caused a modest bias (ΔCBF between -3.78 [P < .001] and 2.83 [P < .001] mL/100g/min). CONCLUSIONS: Arterial spin-labeling imaging is reproducible at both field strengths, and the reproducibility is not significantly correlated with age. Furthermore, 3T tolerates more acquisition parameter variations and allows more extensive optimizations so that 3D and 2D acquisitions can be compared.

Development of a diffusion-weighted MRI protocol for multicentre abdominal imaging and evaluation of the effects of fasting on measurement of apparent diffusion coefficients (ADCs) in healthy liver.

OBJECTIVE: To assess the effect of fasting and eating on estimates of apparent diffusion coefficient (ADC) in the livers of healthy volunteers using a diffusion-weighted MRI protocol with b-values of 100, 500 and 900 s mm(-2) in a multicentre study at 1.5 T. METHODS: 20 volunteers were scanned using 4 clinical 1.5-T MR scanners. Volunteers were scanned after fasting for at least 4 h and after eating a meal; the scans were repeated on a subsequent day. Median ADC estimates were calculated from all pixels in three slices near the centre of the liver. Analysis of variance (ANOVA) was used to assess the difference between ADC estimates in fasted and non-fasted states and between ADC estimates on different days. RESULTS: ANOVA showed no difference between ADC estimates in fasted and non-fasted states (p = 0.8) nor between ADC estimates on different days (p = 0.8). The repeatability of the measurements was good, with coefficients of variation of 5.1% and 4.6% in fasted and non-fasted states, respectively. CONCLUSION: There was no significant difference in ADC estimates between fasted and non-fasted measurements, indicating that the perfusion sensitivity of ADC estimates obtained from b-values of 100, 500 and 900 s mm(-2) is sufficiently low that changes in blood flow in the liver after eating are undetectable beyond the variability in the measurements. ADVANCES IN KNOWLEDGE: Assessment of the effect of prandial state on ADC estimates is critical, in order to determine the appropriate patient preparation for biological validation in clinical trials.

Type 2 diabetes mellitus and obesity in young adults: the extreme phenotype with early cardiovascular dysfunction.

AIM: A pilot study to phenotype young adults (< 40 years) with Type 2 diabetes mellitus. METHODS: Twenty people with Type 2 diabetes (aged 18-40 years), 10 lean and 10 obese control subjects underwent detailed assessment, including tagged cardiac magnetic resonance imaging, inflammatory proteins, lipids, vitamin D and maximal oxygen uptake. Outcomes were compared between the group with Type 2 diabetes and the control group. RESULTS: Mean (standard deviation) age, Type 2 diabetes duration and BMI in the group with Type 2 diabetes were 31.8 (6.6) years, 4.7 (4.0) years and 33.9 (5.8) kg/m(2) respectively. Compared with lean control subjects, those with Type 2 diabetes had more deleterious profiles of hyperlipidaemia, vitamin D deficiency, inflammation and maximal oxygen uptake relative to body mass. However, there was no difference between the group with Type 2 diabetes and the obese control group. The group with Type 2 diabetes had a higher left ventricular mass and a trend towards concentric remodelling compared with the lean control group (P = 0.002, P = 0.052) but not the obese control group (P > 0.05). Peak early diastolic strain rate was reduced in the group with Type 2 diabetes [1.51 (0.24)/s] compared with the lean control [1.97 (0.34)/s, P = 0.001] and obese control [1.78 (0.39)/s, P = 0.042] group. CONCLUSIONS: Young adults with Type 2 diabetes and those with obesity have similar adverse cardiovascular risk profiles, higher left ventricular mass and a trend towards left ventricular concentric remodelling. In addition, those with Type 2 diabetes demonstrate diastolic dysfunction, a known risk marker for future heart failure and mortality.

Multicontrast MR imaging at 7T in multiple sclerosis: highest lesion detection in cortical gray matter with 3D-FLAIR.

BACKGROUND AND PURPOSE: 7T MR imaging has led to improved detection and classification of cortical MS lesions, mainly based on T2*-weighted gradient-echo sequences. Depiction of cortical GM by using the recommended MS imaging protocol has not yet been investigated at 7T. We aimed to investigate prospectively which recommended sequence for clinical use has the highest value at 7T, in terms of GM and WM lesion detection. MATERIALS AND METHODS: Thirty-seven patients with MS (mean age, 43.8 years; 25 women) and 7 healthy controls (mean age, 40.4 years; 5 women) underwent multicontrast 7T MR imaging including the recommended clinical 2D-T2WI, 3D-T1WI, 3D-FLAIR, and GM-specific 3D-DIR. Lesions were scored and categorized anatomically by 3 raters, in consensus. The value of sequences was evaluated lesion-wise and patient-wise (Wilcoxon signed-rank test). RESULTS: At 7T, 3D-FLAIR detected the highest number of total cortical GM lesions (217), 89% more than 3D-DIR and 87% and 224% more than 2D-T2WI and 3D-T1WI. Patient-wise analysis showed that this difference between 3D-FLAIR and 3D-DIR was statistically significant (P<.04), and most pronounced for the number of mixed lesions (P<.03). 3D-FLAIR also detected the highest number of total WM lesions (2605), but the difference with 3D-DIR and 3D-T1WI was not significant. CONCLUSIONS: When using recommended clinical sequences at 7T, the best way to detect cortical GM lesions is with 3D-FLAIR and not by GM-specific 3D-DIR or by conventional 2D-T2WI and 3D-T1WI sequences.

Retinoblastoma: value of dynamic contrast-enhanced MR imaging and correlation with tumor angiogenesis.

BACKGROUND AND PURPOSE: Noninvasive evaluation of retinoblastoma treatment response has become more important due to increased use of eye-sparing treatments. We evaluated the relation between DCE-MR imaging and histopathologic parameters to determine the value of DCE-MR imaging in assessing tumor angiogenesis and prognostic features. MATERIALS AND METHODS: Fifteen consecutive patients with retinoblastoma (mean age, 24 months; range, 2-70 months) undergoing enucleation as the primary treatment (15 eyes) were scanned at 1.5T by using dedicated surface coils. Pretreatment DCE-MR imaging of the most affected eye was evaluated by 2 observers by using curve-pattern analysis, with the first 5 minutes of each curve and the full time-series described as κ(5min) and κ(17min), respectively. Assessed histopathologic and immunologic parameters included optic nerve invasion, choroid invasion, MVD, tumor necrosis, and expression of VEGF and Flt-1. RESULTS: The median value of κ(5min) was 1.28 (range, 0.87-2.07) and correlated positively with MVD (P = .008). The median value of κ(17min) was 1.33 (range, 0.35-3.08) and correlated negatively with tumor necrosis (P = .002). Other histopathologic and immunohistopathologic parameters did not correlate with DCE-MR imaging parameters. Interobserver agreement was 0.53 for κ(5min) and 0.91 for κ(17min). CONCLUSIONS: In retinoblastoma, the early phase of the DCE time curve positively correlates with MVD, while the presence of late enhancement is correlated with necrosis. Thus, the potential for DCE-MR imaging to noninvasively assess tumor angiogenesis and necrosis in retinoblastoma is promising and warrants further investigation.

A study of the brain's resting state based on alpha band power, heart rate and fMRI.

Considering that there are several theoretical reasons why fMRI data is correlated to variations in heart rate, these correlations are explored using experimental resting state data. In particular, the possibility is discussed that the "default network", being a brain area that deactivates during non-specific general tasks, is a hemodynamic effect caused by heart rate variations. Of fifteen healthy controls ECG, EEG and fMRI were co-registered. Slice time dependent heart rate regressors were derived from the ECG data and correlated to fMRI using a linear correlation analysis where the impulse response is estimated from the data. It was found that in most subjects substantial correlations between heart rate variations and fMRI exist, both within the brain and at the ventricles. The brain areas with high correlation to heart rate are different from the "default network" and the response functions deviate from the canonical hemodynamic response function. Furthermore, a general negative correlation was found between heart beat intervals (reverse of heart rate) and alpha power. We interpret this finding by assuming that subject's state varies between drowsiness and wakefulness. Finally, given this large correlation, we re-examined the contribution of heart rate variations to earlier reported fMRI/alpha band correlations, by adding heart rate regressors as confounders. It was found that inclusion of these confounders most often had a negligible effect. From its strong correlation to alpha power, we conclude that the heart rate variations contain important physiological information about subject's resting state. However, it does not provide a full explanation of the behaviour of the "default network". Its application as confounder in fMRI experiments is a relatively small computational effort, but may have a substantial impact in paradigms where heart rate is controlled by the stimulus.

Cerebral blood flow by using pulsed arterial spin-labeling in elderly subjects with white matter hyperintensities.

BACKGROUND AND PURPOSE: On MR imaging, white matter hyperintensities (WMH) on T2-weighted images are generally considered as a surrogate marker of ischemic small vessel disease in elderly subjects. Pulsed arterial spin-labeling (PASL) is a noninvasive MR perfusion-weighted technique. We hypothesized that elderly subjects with diffuse confluent WMH should have lower cerebral blood flow (CBF) measurements than subjects with punctiform or beginning confluent WMH. MATERIALS AND METHODS: MR images of 21 subjects (13 women; mean age, 76 years; SD, 5), stratified for the degree of WMH, from a single center within the multinational Leukoaraiosis and Disability (LADIS) study, were investigated. CBF images were obtained by means of quantitative imaging of perfusion by using a single-subtraction second version, with thin-section TI periodic saturation PASL. Values of cortical gray matter, subcortical (including white matter and deep gray matter), and global CBF were calculated. CBF measurements of subjects with diffuse confluent WMH (n = 7) were compared with those of subjects with punctiform or beginning confluent WMH (n = 14). RESULTS: Subjects with diffuse confluent WMH were found to have approximately 20% lower mean global CBF (43.5 mL/100 mL/min; SD, 6.3) than subjects with punctiform or beginning confluent WMH (57.9 mL/100 mL/min; SD, 8.6; P < .01), as well as approximately 20% lower mean subcortical (P < .01) and cortical gray matter CBF (P < .05). CONCLUSION: PASL revealed a significant reduction of CBF measurements in elderly subjects with diffuse confluent WMH.

Artifact removal in co-registered EEG/fMRI by selective average subtraction.

OBJECTIVE: Co-registration of EEG (electroencephalogram) and fMRI (functional magnetic resonance imaging) remains a challenge due to the large artifacts induced on the EEG by the MR (magnetic resonance) sequence magnetic fields. Thus, we present an algorithm, based on the average-subtraction method, which is able to correct EEG data for gradient and pulse artifacts. METHODS: MR sequence timing parameters are estimated from the EEG data and both slice and volume artifact templates are subtracted from the data. A clustering algorithm is proposed to account for the variability of the pulse artifact. RESULTS: The algorithm is able to keep the spontaneous EEG as well as visual evoked potentials (VEPs), while removing gradient and pulse artifacts with only a subtraction of selectively averaged data. In the frequency domain, the artifact frequencies are strongly attenuated. Estimated MR sequence time parameters showed that the correction is extremely sensitive to the slice time value. Pulse artifact clustering showed that most of the variability is due to the time jitter of the pulse artifact markers. CONCLUSIONS: Selective subtraction of averages in combination with proper time alignment is enough to remove most of the MR-induced artifacts. SIGNIFICANCE: Clean EEG can be obtained from raw signals that are corrupted by MR-induced artifacts during simultaneous EEG-fMRI scanning without using dedicated hardware to synchronize EEG and fMRI clocks.

The hemodynamic response of the alpha rhythm: an EEG/fMRI study.

EEG was recorded during fMRI scanning of 16 normal controls in resting condition with eyes closed. Time variations of the occipital alpha band amplitudes were correlated to the fMRI signal variations to obtain insight into the hemodynamic correlates of the EEG alpha activity. Contrary to earlier studies, no a priori assumptions were made on the expected shape of the alpha band response function (ARF). The ARF of different brain regions and subjects were explored and compared. It was found that: (1) the ARF of the thalamus is mainly positive. (2) The ARFs at the occipital and left and right parietal points are similar in amplitude and timing. (3) The peak time of the thalamus is a few seconds earlier than that of occipital and parietal cortex. (4) No systematic BOLD activity was found preceding the alpha band activity, although in the two subjects with the strongest alpha band power such correlation was present. (5) There is a strong and immediate positive correlation at the eyeball, and a strong negative correlation at the back of the eye. Furthermore, it was found that in one subject the cortical ARF was positive, contrary to the other subjects. Finally, a cluster analysis of the observed ARF, in combination with a Modulated Sine Model (MSM) fit to the estimated ARF, revealed that within the cortex the ARF peak time shows a spatial pattern that may be interpreted as a traveling wave. The spatial pattern of alpha band response function represents the combined effect of local differences in electrical alpha band activity and local differences in the hemodynamic response function (HRF) onto these electrical activities. To disentangle the contributions of both factors, more advanced integration of EEG inverse modeling and hemodynamic response modeling is required in future studies.

Correlating the alpha rhythm to BOLD using simultaneous EEG/fMRI: inter-subject variability.

Simultaneous recording of electroencephalogram/functional magnetic resonance images (EEG/fMRI) was applied to identify blood oxygenation level-dependent (BOLD) changes associated with spontaneous variations of the alpha rhythm, which is considered the hallmark of the brain resting state. The analysis was focused on inter-subject variability associated with the resting state. Data from 7 normal subjects are presented. Confirming earlier findings, three subjects showed a negative correlation between the BOLD signal and the average power time series within the alpha band (8--12 Hz) in extensive areas of the occipital, parietal and frontal lobes. In small thalamic areas, the BOLD signal was positively correlated with the alpha power. For subjects 3 and 4, who displayed two different states during the data acquisition time, it was shown that the corresponding correlation patterns were different, thus demonstrating the state dependency of the results. In subject 5, the changes in BOLD were observed mainly in the frontal and temporal lobes. Subject 6 only showed positive correlations, thus contradicting the negative BOLD alpha power cortical correlations that were found in most subjects. Results suggest that the resting state varies over subjects and, sometimes, even within one subject. As the resting state plays an important role in many fMRI experiments, the inter-subject variability of this state should be addressed when comparing fMRI results from different subjects.

Regional timing of myocardial shortening is related to prestretch from atrial contraction: assessment by high temporal resolution MRI tagging in humans.

Earlier studies have shown substantial nonuniformity in normal left ventricular (LV) myocardial function concerning both the degree of shortening and timing of shortening. We hypothesized that nonuniform LV function may be related to nonuniform prestretch induced by atrial contraction. Eleven healthy human subjects were studied using MRI myocardial tagging and strain analysis. The amount of circumferential prestretch was assessed in 30 LV segments. Prestretch was defined as the difference in strain between end diastole (at ECG R wave) and diastasis. Furthermore, both the degree of shortening (quantified as peak circumferential shortening, peak systolic shortening rate, and amount of postsystolic shortening) and timing of shortening (quantified as the onset time of shortening and time to peak shortening) were assessed. LV prestretch was found to be nonuniform, with the highest values in the lateral wall. The amount of segmental prestretch correlated significantly with peak shortening (r = 0.79), peak shortening rate (r = 0.50), amount of postsystolic shortening (r = 0.67), onset time of shortening (r = -0.57), and time to peak shortening (r = 0.71) (P < 0.001 for each of these relations). These relations may be explained by regional differences in wall stress or by a regional Frank-Starling effect. The correlation between timing of shortening and prestretch demonstrates that mechanical timing is not determined by electrical phenomena alone. In conclusion, regional variation in LV function correlates with the nonuniform prestretch from atrial contraction.

Timing of cardiac contraction in humans mapped by high-temporal-resolution MRI tagging: early onset and late peak of shortening in lateral wall.

Mechanical asynchrony is an important parameter in predicting the response to cardiac resynchronization therapy, but detailed knowledge of cardiac contraction timing in healthy persons is scarce. In this work, timing of cardiac contraction was mapped in 17 healthy subjects with high-temporal-resolution (14 ms) MRI myocardial tagging and strain analysis. Both the onset time of circumferential shortening (T(onset)) in early systole and the time of peak circumferential shortening (T(peak)) at end systole were determined. The onset of shortening width (time needed for 20-90% of the left ventricle to start shortening) was small (35 +/- 9 ms). A distinct spatial pattern for T(onset) was found, with earliest onset in the lateral wall and latest onset in the septum (P = 0.001). Compared with T(onset), T(peak) had a larger width (121 +/- 22 ms) and an opposite spatial pattern, with peak shortening occurring earlier in the septum than in the lateral wall (P < 0.001). Postsystolic shortening (T(peak) later than aortic valve closure; P < 0.05) was observed in 13 of the 30 cardiac segments, mainly in the lateral and basal segments. Shortening in these segments continued 58 +/- 14 ms after aortic valve closure, during which circumferential shortening increased from 16.9 +/- 1.2% to 20.0 +/- 1.5%. Maps of the timing of contraction in normal subjects may serve as a reference in detecting mechanical asynchrony due to intraventricular conduction defects or ischemia.

Interindividual differences of medial temporal lobe activation during encoding in an elderly population studied by fMRI.

Functional MRI (fMRI) is used to study medial temporal lobe (MTL) activation during encoding of new information into memory. In most studies, fMRI data of different subjects are averaged in standard coordinate space. However, interindividual differences in activation can be extensive, reflecting functional heterogeneity. Further, anatomical differences in brain structure cause additional variance and loss of registration accuracy. Such differences in structural and functional MTL characteristics may interfere with the efficiency of averaging data across subjects, and may become more significant with aging and dementia. The current study concerns the analysis of individual differences in MTL activation associated with episodic encoding.Twenty-nine healthy elderly men between 60 and 70 years old performed a simple face encoding task during fMRI scanning. Individual data were analyzed in native space, and compared to the group average in standard space (Talairach and Tournoux).MTL volumes between subjects varied between 6.34 and 11.27 cm(3), and had considerable variation when mapped to standard space. Eighteen of the 29 subjects showed MTL activity and activation patterns varied both in location and size (ranging from 0.11 to 1.78 cm(3)), with the strongest activation in the left posterior part of the MTL. In standard space, no region was significantly activated on a group level at a comparable alpha level. We conclude that while the majority of elderly subjects show MTL activation during episodic encoding of faces, there is considerable structural and functional variability between subjects. Group analysis in standard space may not be appropriate for studies of a complex structure such as the MTL, particularly not in aging and dementia.

Identifying confounds to increase specificity during a "no task condition". Evidence for hippocampal connectivity using fMRI.

Functional MRI can be applied to study connectivity in the brain during a "no task condition." This study focuses on applying a multiple linear regression analysis to identify spurious connectivity caused by confounding factors such as physiologic noise and to separate these from hippocampal connectivity caused by the blood oxygen level dependent (BOLD) signal during a no-task condition. Regressors of interest (hippocampal time courses) as well as regressors of no interest (respiratory signal and cerebrospinal fluid), were included in the analysis, and each yielded a connectivity map. This method was applied at high sampling rate (limited volume, proper physiologic noise sampling), low sampling rate (whole brain scans possible), and at high and low spatial resolution in five healthy control subjects. Regressors of no interest showed specific connectivity patterns, different from hippocampal regressors. The latter showed connectivity between left and right hippocampus. The current study shows successful application of a multiple regression analysis to study connectivity between left and right hippocampus. Both maps of hippocampal connectivity caused by BOLD signal and connectivity caused by spurious signals could be identified.